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Abstract 

A scaling operation on non-negative integers can be defined in a randomised 
way by transforming an integer into the corresponding binomial distribution with 
success probability being the scaling factor. We explore a similar (thinning) op- 
eration defined on counting measures and characterise the corresponding discrete 
stablility property of point processes. It is shown that these processes are exactly 
Cox (doubly stochastic Poisson) processes with strictly stable random intensity 
measures. The paper contains spectral and LePage representations for strictly 
stable measures and characterises some special cases, e.g. independently scattered 
measures. As consequence, spectral representations are provided for the probabil- 
ity generating functional and void probabilities of discrete stable processes. An 
alternative cluster representation for discrete stable processes is also derived us- 
ing the so-called Sibuya point processes that constitute a new family of purely 
random point processes. The obtained results are then applied to explore stable 
random elements in discrete semigroups, where the scaling is defined by means of 
thinning of a point process on the basis of the semigroup. Particular examples 
include discrete stable vectors that generalise the one-dimensional case of discrete 
random variables studied by Steutel and van Harn (1979) and the family of natural 
numbers with the multiplication operation, where the primes form the basis. 

Key words: cluster process, Cox process, discrete semigroup, discrete stability, 
random measure, Sibuya distribution, spectral measure, strict stability, thinning 
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1 Introduction 



Stability for random variables was introduced by Paul Levy and thereafter became one of 
key concepts in Probability. A random vector £ (or its probability law) is called strictly 
a -stable (notation: StaS) if for any positive numbers a and b the following identity is 
satisfied: 

a l/a £' + b l/a £"= {a + b) 1,a £, 

where £ are independent vectors distributed as £ and = denotes equality in dis- 
tribution. Using t = a /(a + b) in the above definition, the stability property can be 
equivalently expressed as 

t i/a £ + (i-t) i/a c = e (i.i) 

for any t G [0, 1]. 

In a recent work [7J the authors of this paper explored the notion of strict stability 
for a much more general situation of random elements taking values in a commuta- 
tive semigroup X. Stability property relies on two basic operations defined on X: the 
semigroup operation: addition and the rescaling which plays a role of multiplication by 
numbers. Many properties of the classical stability still hold for this general situation, 
but there are many notable differences related to specific algebraic properties of X, e.g. 
the validity of distributivity laws, relation between the zero and the neutral elements, 
compactness, etc. Most importantly, the stable random elements arise as a sum of points 
of a Poisson process which intensity measure has a specific scalable form. In the classical 
case of d- dimensional vectors, this result turns into the LePage representation of strictly 
a-stable laws: any such law corresponds to the distribution of the sum of points of a 
Poisson process IT a in M. d with the density (intensity) function having a product form 
9 a (dp)a(ds) in the polar coordinates (p,s). Here a is any finite measure on the unit 
sphere S d_1 called the spectral measure, and a ((r, oo)) = r~ a . The underlying reason is 
that the process U a is itself stable with respect to scaling and superposition, i.e. 

t 1 ^n' a +(i-t) 1 ^ = u a , (1.2) 

this property being granted by the measure 8 a that scales in a 'right' way and the 
superposition property of Poisson processes. This fundamental stability property leads 
to a characterisation of strictly stable laws in very general abelian semigroups forming 
a cone with respect to rescaling with a continuous argument t > 0, see [TJ. 

The case of discrete spaces, however, cannot be treated the same way since the scaling 
by a continuous argument cannot be defined in such spaces. While infinite divisibility 
of random elements in general semigroups, at least in the commutative case, is well 
understood, see [2H E2], a systematic exploration of stable laws on possibly discrete 
semigroups is not available. This prompts to define rescaling as a direct transformation 
of probability distributions rather than being inherited from rescaling of the underlying 
phase space of the random elements. 

Before the current work the family Z + of non-negative integers with addition was 
the only discrete semigroup for which stability was defined. The discrete stability con- 
cept for non-negative integer valued random variables was introduced by Steutel and 
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van Harn [21]. Since multiplication by any real positive t cannot be directly denned on 
Z+, the result of rescaling becomes a probability measure. Namely, the result ton of 
rescaling of n G Z + by t G [0, 1] becomes the Binomial probability measure Bin(n, t) with 
convention that o n = for any n. Recall that Bin(n, t) corresponds to the sum of n 
independent Bernoulli random variables Bin(l, t) with parameter t. Since one can view n 
as the total count of positive integers between 1 and n, a number sampled from Bin(n, t) 
can then be regarded as total count of integers between 1 and n where each number 
is counted with probability t independently of others. In terms of random variables, 
a random variable £ with distribution \x is mapped to a random variable t o £ having 
distribution of Yln=oPm where {f3 n } is a sequence of independent Bernoulli Bin(l,£)- 
distributed random variables. Thus t o £ can be viewed as a doubly-stochastic random 
variable whose distribution depends on realisations of £. 

Alternatively, if g^(s) = E / denotes the probability generating function (p.g.f.) of 
£, then the p.g.f. of t o £ is given by the composition of g% and the p.g.f. of Bernoulli 
Bin(l, t) law: 

gtatis) = gt(l - t(l - s))) . (1.3) 

Similar to (11. ip . a random variable £ is called discrete a-stable if for all < t < 1 
the following identity is satisfied 

tV'o^^-fofij, (1.4) 

where £" are independent copies of £, see [24J. The full characterisation of discrete 
a-stable random variables (denoted by DaS in the sequel) is provided in [21], where it 
is shown that the laws satisfying (II. 4p exist only for a G (0, 1] and each such law has a 
p.g.f. of the form 

Es* = exp{-c(l - s) a } (1.5) 

for some c > 0. 

In this paper we show that the discrete stability of non-negative integer random 
variables is a particular case of stability with respect to the thinning operation of point 
processes defined on a rather general phase space. If the phase space consists of just one 
element, a point process may have multiple points, their total number is a non-negative 
integer random variable and the thinning of the points is equivalent to the above defined 
scaling operation o on discrete random variables. All known properties of discrete stable 
laws have their more general counterparts in point process settings. In particular, we 
show in Section [3] that a-stable point processes with respect to thinning are, in fact, 
Cox processes with the (random) parametric measure being a positive a-stable measure 
with a G (0, 1]. Because of this reason we first address general strictly stable random 
measures in Section [21 Note in this relation that so far only independently scattered 
stable measures received much attention in the literature, see [23] and more recently [B] 
and [12]. 

The central result is Theorem 13 . 31 which shows that discrete stable point processes are 
exactly Cox processes driven by a strictly stable random measure on the phase space. It 
should be noted that Cox processes driven by various random measures are often used in 
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spatial statistics, see [US, EH1 EE] . A particular feature of discrete stable processes is that 
the point counts have discrete a-stable distributions, which have infinite expectations 
unless a = 1. 

It is known that discrete stable random variables can be represented as the sum of 
Poisson number of Sibuya distributed integer random variables, see [9]. In Section 0] we 
show that an analogous cluster representation of finite discrete stable point processes 
also holds. The clusters are Sibuya point processes which seem to be a new class of point 
processes not considered so far. 

Some of our results for general point processes, especially defined on a non-compact 
spaces do not have, however, counterparts for discrete random variables or become triv- 
ial. This concerns the results of Section [5] which draws analogy with infinitely divisible 
point processes. 

Another important model arising from the general point process setting is the case 
of random elements in discrete semigroups which possess at most countable basis. We 
show in Section [6] that the uniqueness of representation of each element as a linear finite 
combination with natural coefficients of the basis elements is a necessary condition to 
be able to define discrete stability. We give characterisation of DaS elements in these 
semigroups and establish discrete analogue of their LePage representation. 

The presentation of the theoretical material is complemented by examples from ran- 
dom measures, point processes and discrete semigroups all over the text. In particular, 
in Section Owe introduce a new concept of multiplicatively stable natural numbers which 
could be interesting on its own and characterise the corresponding distributions. 

2 Strictly stable random measures 

Let X be a locally compact second countable space with the Borel cx-algebra B. The 
family of Radon measures on B is denoted by M. These are measures which have finite 
values on the family Bo of relatively compact Borel sets, i.e. the sets with compact 
topological closure. Note that all measures considered in this paper are assumed to be 
non- negative. The zero measure is denoted by 0. 

A random measure is a random element in the measurable space [M, Ai], where the 
cr-algebra M. is generated by the following system of sets: 



The distribution of a random measure ( can be characterised by its finite-dimensional 
distributions, i.e. the distributions of (((Bi), . . . , ((B n )) for each finite collection of 
disjoint sets Bi, ... , B n G Bo. 

It is well known [51 p. 192] that the distribution of ( is uniquely determined by its 
Laplace functional 



defined on non-negative bounded measurable functions h : X i— > M + with compact 



{/ieU: fi{Bi)<ti, i = l,...,n}, B t e B, u > 0. 
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support (denoted h G BM(Jf)). From now on we adapt a shorter notation 

(h,/j,)= / h(x) fi(dx) , /iGM. 

The family M can be endowed with the operation of addition and multiplication by 
non-negative numbers as 

(/ii + /i 2 )(-) = / i 1 (-) + /i 2 (-); (2.2) 

(^)O =*//(•), t>0. (2.3) 

If ^ is a finite measure, it is possible to normalise it by dividing it by its total mass 
and so arriving at a probability measure. The normalisation procedure can be extended 
to all locally finite measures as follows. Let Bi,B 2 ,... be a fixed countable base of the 
topology on X that consists of relatively compact sets. Append -Bo = X to this base. 
For each /iGM consider the sequence of its values (fi(B ) , fi(Bi) , fi(B 2 ) , . . .) possibly 
starting with infinity, but otherwise finite. Let be the smallest non-negative integer 
i for which < [i{Bi) < oo, e.g., i(fi) = if fi is a finite measure. The set 

§ = {// G M : n(B i{jl) ) = 1} (2.4) 

is measurable, since 

oo 

§ = Mi U |J {fi G M : /i(B ) = oo, M^x) = • • • = M^n-i) = 0, /i(B n ) = 1} , 

n=l 

where Mi is the family of probability measures on X. Note that Sfl {/i : fi(X) < oo} = 
Mi. Furthermore, every fi can be uniquely associated with the pair (p,, ^(B^)) G 
§ x R + , so that n = ^{Bi^jl and M = S x R + . 

Definition 2.1. A random measure ( is called strictly a-stable (notation StaS) if 

t 1/a C' + (l-t) 1/a C" = C (2- 5 ) 
for all < t < 1, where (" are independent copies of the random measure (. 

Definition 12.11 yields that for any B\, . . . , B n G Bo, the vector (£(£?i), . . . , ((B n )) is 
a non-negative StctS n-dimensional random vector implying that a G (0,1]. It is well 
known that one-sided strictly stable laws corresponding to a = 1 are degenerated, so 
that StaS measures with a = 1 are deterministic. 



Theorem 2.2. A locally finite random measure ( is StaS if and only if ( is deterministic 
in the case a = 1 and in the case a G (0, 1) if and only if its Laplace functional is given 
by 

L c [h) =exp{- / (1 - e- M )A(dfx)\ , h G BM(Jf) , (2.6) 

iM\{0} J 
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where A is a Levy measure, i.e. a Radon measure on M \ {0} such that 

f (l- e - M )A(dfi) < oo (2.7) 

JM\{0} 

for all h G BM(X), and A is homogeneous, i.e. A(tA) = t~ a A(A) for all measurable 
A C M\ {0} and t > 0. 

Proof. Sufficiency. It is obvious that each deterministic measure is StaS with a = 1. We 
show that (12.61) defines a Laplace functional of a random measure. Consider hi, . . . , hj. G 
BM(X) and the image A of A under the map \i i— > ((hi, //),..., (h k , a)). It is easy to see 
that A is a homogeneous measure on Then ^E i=1 ij^i] as a function of ti, . . . , 
is the Laplace transform of a totally skewed to the right (i.e. almost surely positive) 
strictly stable random vector ((hi, (),... , (hk, 0) with the Levy measure A. 

Furthermore, L^[0] = 1 and L^[h n ] — ► if /i n T ^ pointwise. To see this, H n (a) = 
(h n , a) | (h, u) = H(u) by the monotone convergence theorem. In its turn, (H n ,A) 
converges to (H, A) again by the monotone convergence. Now by an analogue of [5J 
Th. 6. 4. II] for the Laplace functional, (I2.6P indeed is the Laplace functional of a random 
measure. 

Necessity. Strictly stable random measures can be treated by means of the general 
theory of strictly stable laws on semigroups developed in [7J. Consider the cone M of 
locally finite measures with the addition and scaling operation defined by (12.21) and (12.31) . 
In terms of [7] this set becomes a pointed cone with the origin and neutral element being 
the zero measure, where also the second distributivity law holds. The gauge function 
\\fi\\ = fi(Bi^), fj, e M, is homogeneous with respect to multiplication of measures by a 
number, so that M = §x M + becomes the polar decomposition of M. As in Ex. 8.6], 
we argue that the characters on M 

xM=exp{-(h,ix)}, heBM(X), 

have a strictly separating countable subfamily and are continuous, so that condition (C) 
of [7] holds. Furthermore, Xh(s^) — > 1 as s j for all fi G M, so that condition (E) 
of [7j is also satisfied. By [?l Th. 5.16], the log-Laplace transform of a StaS measure 
satisfies L ( [sh] = s a L ( [h] for s > 0, where a G (0, 1] by [7, Th. 5.20(h)]. If a = 1, this 
identity implies that the values of ( are deterministic. By [3 Th. 6.5(h), Th. 6.7(i)], the 
log-Laplace functional of a StaS random measure with a G (0, 1) can be represented as 
the integral similar to (12. 6p with the integration taken with respect to the Levy measure 
of ( over the second dual semigroup, i.e. the family of all characters acting on functions 
from BM(X). 

In order to reduce the integration domain to M \ {0} it suffices to show that the 
convergence (h, /i^) — > g(h) for all h G BM(X) implies that g(h) = (h,fi) for some 
H G M. Setting h — Ib above implies convergence and hence boundedness of the 
sequence {fi k (B), k > 1} for each B G B . By [51 Cor. A2.6.V] the sequence {/z^, k > 
1} is relatively compact in the vague topology. Given that (h, Hk) converges for each 
h G BM(X), we obtain that all convergent subsequences of {/ifc, k > 1} share the same 
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limit that can be denoted by fi and used to represent g(h) = (h,fi). Indeed, if two 
subsequences /x^, /ij 2 , . . . and /i^, fMj 2 , . . . have different limits, say, /i' and /i", then there 
is h G BM(X) such that (h, //) ^ (h,fi"). Then for the subsequence /i^ , /i^ , /ij 2 , /Xj 2 , . . . 
(possibly after removing the repeating members) the limit of the integrals of h does not 
exist contradicting the assumption. Now (I2.6P follows from Th. 7.7]. □ 

The properties of StaS measure ( are determined by its Levy measure A. For in- 
stance, if X = W 1 , then ( is stationary if and only if A is invariant with respect to the 
map yu(-) 1 — > fi(- + a) for all a G R d and all /i from the support of A. Note also that (12. 7p 
is equivalent to 

f (x - e -^ B ))A(dfi) < 00 (2.8) 

iM\{0} 

for all B G B . 

It is also useful to decompose A into the radial and directional components using 
the polar decomposition of M described above. Note that the map /1 h-* (ft, fi(Bi^)) is 
measurable. If 

a (A) = A({tfi : \i G A, t > 1}) 

for all measurable Ac§, then A(A x [a, b]) = a(A)(a~ a — b~ a ), so that A is represented 
as the product of a and the radial component given by the measure 8 a defined as 
9 a {[r, 00)) = r~ a , r > 0. By the reason which will become apparent in the proof of 
Theorem 12.31 below, it is more convenient to scale a by the value T(l — a) of the Gamma 
function. The measure a = T(l — a)a is called the spectral measure of (. Note that o 
is not necessarily finite unless X is compact. 

Condition (12.81) can be reformulated for the spectral measure a as 

/ fi{B) a a{d^) < 00 (2.9) 

for all B E Bq- Note that a can also be defined on any other reference sphere §' C M, 
provided each fi G M can be uniquely represented as tfi' for some /i' G S'. 

Theorem 2.3. Let a be the spectral measure of a StaS random measure ( with Levy 
measure A and let a G (0, 1). Then 

L c [h] = exp{- J \h, /i)V(cZ^)} , heBM(X). (2.10) 

Furthermore, 

(i) ( is a.s. finite if and only if its Levy measure A (resp. spectral measure a) is 
supported by finite measures and cr(§) = cr (Mi) is finite; 

(ii) the Laplace functional \2. 6)) defines a non-random measure if and only if a = 1. 
In this case ( = Ti(-) = J s /i(-)cr(d/x). 
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Proof. The representation f 1 2 . 1 1) follows from (12.61) 



since 



(1 _ e -M)A(dn) = / / (1-e 



tM )6 a (dt)a(dfi) 




r 



(1 - a) / (h,fi) a a(dfi) 



Js 



implying (I2.10p . 

(i) Taking h n = lx n , ti> 1, where X n G B form a nested sequence of relatively compact 
sets such that X = U n X n , we obtain that the Laplace transform of ((X) is 



where the integral is finite if and only if A is supported by finite measures. If this is the 
case, the spectral measure a is defined on Mi and 



defines a finite random variable if and only if a (Mi) is finite. 

(ii) If ( is StaS with a — 1, then its values on all relatively compact sets are deterministic, 
and so ( is a deterministic measure. Furthermore, (12.51) clearly holds with a = 1 for any 
deterministic (. Finally, fl 2 . 1 [) with a = 1 yields 



If the spectral measure is degenerate, the corresponding StaS measure has a partic- 
ularly simple structure. 

Theorem 2.4. Assume that the spectral measure a of a StaS measure £ with < a < 1 
is concentrated on a single measure so that a = cS^ for some [i G S and c > 0. Then 
( can be represented as ( = c 1//a C a /i, where ( a is a positive StaS random variable with 
Laplace transform Ee" 2 ^ a = exp{— z a }. 

Proof. It suffices to verify that the Laplace transforms of both measures coincide and 
equal L^[h] = exp{—c(h,fi) a }. □ 

Example 2.5. Let X = IR d and let a be concentrated on the Lebesgue measure I. By 
Theorem 12.41 the corresponding StaS measure ( is stationary and proportional to ( a £. 

If X is compact, then all measures [i G M are finite and condition ( 12.9f) implies 
that <t(Mi) is finite, so that the spectral measure is defined on the family of probability 
measures on X. 




M\{0} 



limL c [/i n ] = exp{-(j(Mi)} 




□ 
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Example 2.6 (Finite phase space). In the special case of a finite X = {x±, . . . ,Xd} the 
family of probability measures becomes the unit simplex = {x G : J^Xj = 1} 
and a StaS random measure ( is nothing else but a <i-dimensional totally skewed (or 
one-sided) strictly stable random vector C = (&>•••> Cn)- Its Laplace transform is given 
by 

Ee- {h '° = ex p{~y (M)M^)J , (2.11) 

where a is a finite measure on the unit simplex. Alternatively, the integration can be 
taken over the unit sphere. It is shown in [T7j that this Laplace functional can be written 

as 

Ee -<ft,0 = e -H K (h«) ; heR d + , (2.12) 

where h a = (hf, . . . , h%) and Hk is the support function of a certain convex set K that 
appears to be a generalisation of zonoids. 

Example 2.7 (StaS random measures on M. d ). Fix a Radon measure fi on X = M. d and 
let a be the image of a locally finite measure v on M. d under the map x i— > + x). The 
corresponding StaS random measure £ has the Laplace transform 

L ( [h] = exp{- / ( / h(x - y)/j,(dx)) u(dy)\ . 

To ensure condition (12.91) it is necessarily to assume that 

/ (/x(S + y))«i/(dy)<oo (2.13) 

for all relatively compact B. In particular, 

Ee- Z((B) =exp{-z a J (ji(B + ar))° , 5 G i3 . (2.14) 

If z/ is the Lebesgue measure, this construction can be used to define a rather rich family 
of stationary StaS random measures. 

Note that usually in the literature, in particular, in [23| Sec. 3.3], the term a-stable 
measure is reserved for an independently scattered measure, i.e. a measure with inde- 
pendent a-stable values on disjoint sets. Our notion is most general as the following 
proposition shows. 

Theorem 2.8. A StaS random measure with a G (0, 1) is independently scattered if 
and only if its spectral measure a is supported by the set {5 X : x G X} C § of Dirac 
measures. 

Proof. Sufficiency. The Laplace functional ( 12. 10D of ( becomes 

L c [h] = exp j- I h a {x) a(dx)} , 
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where a is the image of a under the map S x i— > x. It is easy to see that i<[X)iLi hi] = 
IliLi-^cM f° r a U functions hi,...,h n with disjoint supports, meaning, in particular, 
that for any disjoint Borel sets Bi, ... , i? n the variables C(-^i) with Laplace transforms 
L(\zi IbJ, i = 1, . . . ,n, are independent, so that ( is independently scattered. 
Necessity. Since C(-^i) an d C(B 2 ) are independent for disjoint Bi, B 2 G £>, we have that 

L c [l Bl + l B2 } = L c [l Bl }L c [l B2 }. 

By (jno]), 

/ ^(Sx) + //*(5 2 ) - (/i(-Bi) + /i(5 2 )) Q ] a(^) = . 

Since a G (0, 1), the integrand is a non-negative expression whatever //(-Bi), fJ-(B 2 ) > 
are. Since a is a positive measure, the integral is zero only if the integrand vanishes 
on the support of a, i.e. for a-almost all \x and all disjoint Bi,B 2 either fi(Bi) = or 
n(B 2 ) = meaning that \i is concentrated at a single point. □ 

Along the same lines it is possible to prove the following result. 

Theorem 2.9. The values C{Bi), ■ ■ ■ , C{B n ) of a StaS random measure ( with a G (0, 1) 
on disjoint sets Bi, . . . , B n are independent if and only if the support of the spectral 
measure a (or of the Levy measure A ) does not include any measure that has positive 
values on at least two sets from B ly . . . , B n . 

The following theorem provides a LePage representation of a StaS measure. 
Theorem 2.10. A random measure ( is StaS if and only if 

C=X>' (2-15) 

we* 

where is a Poisson process on M\ {0} driven by a measure A satisfying A2.7\) and such 
that A(tA) = t~ a A(A) for allt>0 and any measurable A. In this case A is exactly the 
Levy measure of (. 

If the spectral measure a corresponding to A satisfies c = < oo, then 

oo 

( = 65>;". S4 , h= (^—\\ (2.16) 

k=l v ; 

where S\, e 2 , . . . are i.i.d. random measures with distribution c~ l a and 7^ = £1 + ■ • ■ + 
k > 1, for a sequence of independent exponentially distributed random variables with 
mean one. 

Proof. Consider a continuous function h from BM(A) and define a map u. 1— > (h, u) 
from M to M + = [0, 00). The image of ^ under this map becomes a Poisson process 
{xi, i > 1} on M + with intensity measure 9 that satisfies 

0([r,oo))=A(Oi: (h,n)>r}). 
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Since l—e~ x > al x > r for any r > and some constant a > 0, Condition ( 12. 7ft implies that 
the right-hand side is finite. The homogeneity property of A yields that 9([r, oo)) = cr~ a 
for all r > and some constant c > 0. 

It is well known (see [23] and [7]) that for such intensity measure 9 the sum ^ Xj 
converges almost surely. Thus, J] e9 (/i, /ij) converges for each continuous ft, G BM(Jf), 
so that the series (12.151) converges in the vague topology to a random measure, see J2UJ 
Eq. (3.14)]. 

By [20l Prop. 3.19] a sequence of random measures converges weakly if and only if 
the values of their Laplace functionals on any continuous function converge. This is seen 
by noticing that the Laplace functional ( 12.61) of ( coincides with the Laplace functional 
of the right-hand side of (12. 15ft . where the latter can be computed as the p.g.fl. of a 
Poisson process, see (13. 2ft . 

Finally, (I2.16P follows from the polar decomposition for A. □ 

3 Discrete stability for point processes 

Point processes are random elements with realisations in the set of locally finite counting 
measures. Each counting measure (p can be represented as the sum (p = ^2S Xi of 
unit masses, where we allow for the multiplicity of the support points {xj, % > 1}. 
The distribution of a point process $ can be characterised by its probability generating 
functional (p.g.fl.) defined on functions u : X h (0, 1] such that 1 — u G BM(Jf) by 
means of 

G 9 [u] = Eexp{(logw,$)} = E Yl u(xif i{x * }) = L#[-logu] (3.1) 

with convention that the product is 1 if $ is a null-measure (i.e. the realisation of the 
process contains no points). The p.g.fl. can be extended to pointwise monotone limits 
of the functions from BM(X) at the expense of assuming infinite and zero values. 

A Poisson process with intensity measure \x is characterised by the property that 
for any finite collection of disjoint sets Bi, . . . , B n G Bq the variables U^Bi), . . . , n M (S„) 
are mutually independents Poisson distributed random variables with means /j,(Bi), . . . , fi(B n ). 
The p.g.fl. of a Poisson process I1 M is given by 

G n >] = exp{-<l-u,/z)}. (3.2) 

If \x itself is a random measure (, then the obtained (doubly stochastic) point process is 
called a Cox process. Its p.g.fl. is given by 

Gn ( [u} = Eex V {-(l-u,C)}=L ( [l-u]. (3.3) 

Note that the Cox process is stationary if and only if the random measure ( is stationary. 

Addition of counting measures is well defined and leads to the definition of the 
superposition operation for point processes. However, multiplication of counting measure 
by positive numbers is not possible to define by multiplying its values — they no longer 
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remain integers for arbitrary multiplication factors. In what follows we define a stochastic 
multiplication operation that corresponds to the thinning operation for point processes. 
Namely each unit mass 5 Xi in the representation of the counting measure (p = Yli 5 Xi is 
removed with probability 1 — t and retained with probability t independently of other 
masses. The resulting counting measure toip becomes random (even if (p is deterministic) 
and is known under the name of independent thinning in the point process literature, 
see, e.g., [5J, Ch. 8.2] or [161 C. 7], where relations to cluster and Cox processes were 
established. 

To complement this pathwise description, it is possible to define the thinning oper- 
ation on probability distributions of point processes. Namely the thinned process to$ 
has the p.g.fl. 

G t0 $[u] = G*[tu + 1 - 1] = G*[l - t(l - u)} , (3.4) 

see, e.g., (5, p. 240]. From this, it is easy to establish the following properties of the 
thinning operation. 

Theorem 3.1. The thinning operation o is associative, commutative and distributive 
with respect to superposition of point processes, i.e. 

ti o (t 2 o $) = (tit 2 ) o $ = t 2 o (ti o $) , 
t o ($ + $') = (t o $) + (t o $') . 

for any t, t\,t 2 G [0, 1] and any independent point processes $ and 
For disjoint Borel sets B\ and B 2 random variables (t o $)(5i) and (t o $)(5 2 ) are 
conditionally independent given $ and there exists a coupling of $ and t o $ (described 
above) such that (t o $)(-B) < $(-B) almost surely for any B £ B. 

If the phase space X consists of one point, the point process $ becomes a non-negative 
integer random variable <&{X). Since (13.41) turns into fll.3|) . the thinning operation 
becomes the classical discrete multiplication operation acting on positive integer random 
variables. Keeping this in mind, we can define the notion of discrete stability for point 
processes. 

Definition 3.2. A point process $ (or its probability distribution) is called discrete 
a-stable or a-stable with respect to thinning (notation DaS), if for any < t < 1 one 
has 

i 1 /* o $' + (i _ t)V« o £ $ , (3.5) 
where $' and $" are independent copies of $. 

Definition 13.21 and (13.41) yield that $ is DaS if and only if its p.g.fl. possesses the 
property 

G$[l-t 1/Q (l-n)]G'$[l-(l-t) 1/a (l-u)] =G*[u], (3.6) 

for any positive function it such that 1-uG BM(X). 

Relation (13. 5p implies that $(5) is a discrete a-stable random variable for any B G B. 
Hence DaS point processes exist only for a G (0, 1] due to already mentioned result (11.51) . 
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Theorem 3.3. A point process $ is DaS if and only if it is a Cox process IT^ with StaS 
intensity measure (. 

Proof. Sufficiency. If ( is StaS random measure, then the corresponding Cox process 
has the p.g.fi. 

Gu ( [u] = L ( [l —u] = exp{- ! (1 - e- {1 - u ^)A{dfi)\ . (3.7) 

Jm\{o} J 

Thus, 

- t{l - u)] = exp{- / (1 - e- {1 - u ' t ^ ) )A(dfx)} 

Jm\{o} J 

= exp{-t Q / (l-e-^-^^A^/i)) 

by the homogeneity of A. Then (13. 6p clearly holds. 
Necessity. By iterating (13. 5p m times we arrive at 

m' l/a o $! + ••• + m" 1/a o $ m = $ 

for i.i.d. $i,...,$ m implying that DaS point processes are necessarily infinitely 
divisible and 

(G$[l -m" 1/Q (l -w)]) m = G 9 [u] (3.8) 

for all m > 2. 

The crucial step of the proof aims to show that the functional 

L[u] = G 9 [l-u] , mGBM(X), (3.9) 

is the Laplace functional of a StaS random measure (. While the functional L in the 
left-hand side should be defined on all (bounded) functions with compact supports, it is 
apparent that the p.g.fi. G$ in the right-hand side of (13.91) may not be well defined on 
1 — u. Indeed, 1 — u does not necessarily take values from the unit interval. To overcome 
this difficulty, we employ (13.81) and define 

L[u] = - m~ 1/a u}) m . (3.10) 

Since u E BM(X), for sufficiently large m the function 1 — vrT x l a u takes values in [0, 1] 
and equals 1 outside a compact set. Since (I3.10P holds for all m, it is possible to pass 
to the limit as m -» oo to see that 

L[u] = exp < — lim m(l — G$[l — m~ l l a u}) \ . 

L m—*oo J 

By the Schoenberg theorem [31 Th. 3.2.2], L is positive definite if 1 — G$[l — m~ l l a u] is 
negative definite, i.e. 

n 

0^(1 - C # [l - m- 1/a ( Ui + Uj)}) < 
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for all n > 2, u±, . . . , u n G BM(X) and C\, . . . , c n with q = 0. In view of the latter 
condition, we need to show the inequality 

n 

CiCjG*[l - m- 1/a { Ui + Uj)} > . 



Note that 



1 — m l / a u 
lim — = 1 

m — >oo C~ m u 



Denote v j = e~ m u \ i — 1, . . . , n. The continuity of the p.g.fl. G$ now implies that it 
suffices to check that 

n 

*,i=i 

which is exactly the positive definiteness of Gq>. 

Thus, i([Ei=i^»] as a function of ti, . . . , > is the Laplace transform of a 
random vector. Arguing as in the proof of sufficiency in Theorem I2.2[ it is easy to check 
the continuity of L, so that L is indeed the Laplace functional of a random measure (. 
Condition (13.61) rewritten for L means that ( is StaS. □ 

By Theorem 12.31 a StaS random measure ( has the characteristic exponent a = 1 
if and only if £ is deterministic. Respectively, a DaS processes with a = 1 is a Poisson 
process driven by (non-random) intensity measure (. 

Using decomposition A = a x 9 a , the first identity in (13 .7p and (I2.10p . we obtain the 
following result. 

Corollary 3.4 (Spectral representation). A point process $ is DaS with a G (0,1] if 
and only if its p.g.fl. has the form 

G^[u] = expj- J{1 - w,/i)V(d/i)| , I-mgBM(X) (3.11) 

for some locally finite spectral measure a on § that satisfies \2. 

The number of points $(-B) of a DaS process $ in a relatively compact Borel set B 
has the p.g.f. 

Es $(B) =exp|-(l-s) a ^/i(5)V(d/i)} (3.12) 

and so $(5) is either zero or has infinite expectation for < a < 1, while E$(S) is 
finite in the Poisson case a = 1. Furthermore, the avoidance probability is given by 



P{$(£) = 0} = expj- J^{B) a 



so that (12.91) guarantees that the avoidance probabilities are positive. 
Theorem 13.31 immediately implies the following result. 
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Corollary 3.5 (LePage representation of a DaS process). A DaS process $ with Levy 
measure A can be represented as a superposition of (dependent) Cox processes: 

$=E n ^ (3.13) 

where $ is a Poisson process on M \ {0} with intensity measure A. In particular, if its 
spectral measure a is finite then 

oo 

$ = E n , 7fe -^' ( 3 - 14 ) 

k=l 

where Sk, Jk and b are defined in Theorem 12.101 and the sum in A3.14\) is almost surely 
finite. 

Proof. By Theorem 13. 3[ $ is a Cox process IT^, where £ is representable as (12.151) . 
Conditioned on a realisation of we have that IT^ = ^ 11^. by the superposition 
property of a Poisson process. Now the statement easily follows. 

The finiteness of the sum in (13. 14ft follows from the Borel-Cantelli lemma, since 

P^-v^PO > I 7fc } = 1 -exp{-6 7 - 1/a } = 0{k-^) 

for almost all realisations of {7^}. □ 

Corollary 3.6 (see [TO] ) . A random variable £ with non-negative integer values is dis- 
crete a-stable if and only if £ is distributed according to the Poisson law which mean is 
given by a positive strictly stable random variable. The p.g.f of £ is given by 

E^ = e- c(1 - s)a , «e (0,1], 

for some c > 0. 

Example 3.7 (Discrete stable vectors, see Example 12.61) . A point process $ on a finite 
space X can be described by a random vector (£1, •••,£<*) of dimension d with non- 
negative integer components. If $ is DaS then (£i,...,£d) are said to form a DaS 
random vector. For instance, if $ is a DaS point process in any space X, then the point 
counts <&(f?i), . . . , $(-B n ) form a DaS random vector for B 1 , . . . , B n e £> . 

Theorem 13.31 implies that £ = (£1, ...,£<*) is DaS if and only if its components are 
Poisson distributed with means £ = (d, • • • , Cd) an d are conditionally independent given 
£, where £ is a strictly stable non- negative random vector with the Laplace transform 
( 12. lip . Notice that in general the components of £ are dependent DaS random variables, 
unless the spectral measure a of £ is supported by the vertices of the simplex only. 
In view of ( 12. 12ft the p.g.f. of £ can be represented as 

E]Jsf = exp{- f (l-8,^) a a(dfA=e- HK ^-'^, 
i=i I J^a J 

where s = (s\, . . . ,Sd) and 1 = (1, . . . , 1). Thus, the distribution of £ is determined by 
the values of the support function of Hk and its derivatives at 1. 
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Example 3.8 (Mixed Poisson process with stable intensity, see Example 12 .5D . Let X = R d 
and let a attach a mass c to the measure ai on M. d , where I is the Lebesgue measure 
and a > is chosen so that ai G §. Then G<s>[u] = exp{— ca a (l — u,£) a } which is the 
p.g.fl. of a stationary Cox process with a StaS density, i.e. a stationary Poisson process 
in W 1 driven by random intensity measure ac 1 ^ a ( a £, where ( a is a positive StaS random 
variable with Laplace transform E e~ zC - a = e~ z " . This type of Cox process is also known 
as a mixed Poisson process. In particular, if c = a _1// °, then Es _$ ' B ' = e~ ( - 1 ~ s ^ c ' i ^ a , so 
that P{$(Z?) = 0} = e~^ B ^ a for all B G Bq. Note that this point process has infinite 
intensity measure. 

The following construction provides a rather general family of DaS processes. 

Example 3.9 (see Example 12.71) . Consider StaS random measure ( with the spectral 
measure being the image of a locally finite measure v on IR d under the map x t— > fi(- + x) 
for a given Radon measure fi satisfying f 1 2 . 1 3 j) . The corresponding DaS process $ has 
the p.g.fl. 

G$[n] = exp{- / ( / (1 - u(x - y))fi(dx)) u(dy)\ 
and the avoidance probabilities for $ are given by 

P{$( J B) = 0} = expj-y (ji(B + x)Y V(da;)} 

If v is the Lebesgue measure, this construction defines a family of stationary DaS pro- 
cesses parametrised by a measure \i. 

Discrete a-stable point processes appear naturally as limits for thinned superpositions 

of point processes. Let \l/ be a point process on X and let S n — H h ^ n be the sum 

of independent copies of The p.g.fl. of the thinned point process a n o S n is given by 

G 0n o5„[l - /i] = - a n /i]) n , 

where a„ — > is a certain normalisation. On the other hand, the superposition can 
be normalised by scaling its values as a n S n , so that the result of scaling is a random 
measure, but no longer counting. The following basic result establishes a correspondence 
between the convergence of the thinned and the scaled superpositions. 

Theorem 3.10. Let S n , n > 1, be point processes. Then for some sequence {a n }, 
a n S n weakly converges to a non-trivial random measure which is necessarily StaS with a 
spectral measure a if and only if a n o S n weakly converges to a non-trivial point process 
which is necessarily DaS with a spectral measure a. 

The measure a is defined via its finite-dimensional distributions: if £ = (^(f^), . . . , ^/(B d )) 
for Bi, . . . , B d G B , then 

(^ 1 ,..., Md ))eA})=r(l-a)limnP{e/||e||ieA; ||£||i > b n } (3.15) 

n— >oo 

for all measurable A from the unit £i-sphere {x G IRl : ||rc||i = 1}. 
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Proof. The equivalence of the convergence statements is established in [T4"l Theorem 8.4], 
see also (5} Th. 9.3.III]. 

Finally, (13 . 15[) is the standard condition for a n S n to have a limit which is valid in much 
more general setting than for random measures, see Th.4.3] and [Ij. The Gamma- 
term stems from the normalisation of the spectral measure adopted in Section [2J □ 

If X consists of one point, then we recover the result of jl] concerning a relationship 
between discrete stable random variables and conventional stable laws. If a point process 
\1/ has a finite intensity measure, then the strong law of large numbers implies that the 
limit of n~~ 1 S n is deterministic, so that the limiting DaS process is, in fact, Poisson. 
Therefore a DaS limit with a G (0, 1) is only possible for a point process ^ with infinite 
intensity. 

Remark 3.11. Note that not all point processes can give a non-trivial limit in the above 
schemes. Consider, for instance, a point process \1/ defined on J = {1,2,...} with the 
point multiplicities being independent discrete stable with characteristic exponents 
oij = 1/2 + l/(2i), i = 1, 2, . . . . Then n~ l l a o S n for a < 1/2 gives null-process as the 
limit, while for larger a the limiting process is infinite on all sufficiently large i. 

4 Sibuya point process and cluster representation 

Recall that a general cluster process is defined by means of a centre point process N c 
in some phase space Y and a countable family of independent daughter point processes 
N(-\y) in a phase space X indexed by the points of Y. Their superposition in X defines 
a cluster process. The p.g.fl. G[h] of a cluster process is then a composition G c [G4/i|-]] 
of the centre and the daughter processes, see, e.g., p3, Prop. 8. 2. II]. 

Noting (13.21) . the form of p.g.fl. (13.71) suggests that a DaS process can be regarded 
as a cluster process with centre processes being Poisson with intensity measure A in 
Y = M \ {0} and daughter processes being also Poisson parametrised by their intensity 
measure fi G suppA. We shall embark on exploration of the general case in the next 
section, but here we concentrate on the case when measure A charges only finite measures 
and give alternative explicit cluster characterisation of such DaS processes. 

Definition 4.1. Let fi be a probability measure on X. A point process T on X defined 
by the p.g.fl. 

G r [u] = GrwM = 1 - (1 -u,fx) a (4.1) 

is called the Sibuya point process with exponent a and parameter measure jj,. Its distri- 
bution is denoted by Sib(a,/i). 

If X consists of one point, then the point multiplicity is a random variable rj with 
the probability generating functional 

E^ = 1-(1-^) Q , z 6(0,1]. 

In this case we say that rj has the Sibuya distribution denoted by Sib(a). The Sibuya 
distribution corresponds to the number of trials to get the first success in a series of 



17 



Bernoulli trials with probability of success in the kth trial being a/k, see also [9] and pU] 
for efficient algorithms of its simulation. 

In particular, (14.11) implies that T(B) with B G Bq has p.g.f. 

Eu r{B) = 1 -fi a (B)(l -u) a . 

Note that T(B) is non-integrable if fi{B) does not vanish. Furthermore, 

P{T(fl) = 0} = l-/i%B) s 

P{T(B) = 1} = «/i Q (5) = qi {a) //*(£) , 

P{T(X) = n} = (1 - a) (1 - £) . . . (1 _ _2L.) 2 d 4 f fe(a) ^ n > 2 

2 ra — 1 n 

Therefore, conditioned to be non-zero, T(JB) has Sibuya distribution with parameter a 
justifying the chosen name for this process T. In particular T(X) is non-zero a.s. and 
follows the Sibuya distribution. 

Developing the p.g.fl. (14.11) makes it possible to get insight into the structure of a 
Sibuya process: 

oo 
n=l 

u(x n ) (J,(dxi) ■ ■ ■ fi(dx n ) . 

Therefore, as we have already seen, the total number of points of T follows Sibuya distri- 
bution and given this total number, the points are independently identically distributed 
according to the distribution \i. This also justifies the fact that (14. ip indeed is a p.g.fl. 
of a point process constructed this way. This type of processes is called purely random 
in [TBI P- 104]. 

Assume now that the Levy measure A of a StaS random measure ( is supported by 
finite measures. Then (12.101) holds with the spectral measure a defined on Mi so that 
the p.g.fl. of the corresponding DaS process takes the form 

C$N = exp{ / (G T (m)M " iMdfi)} (4.2) 

with Gr(fj,)[u] given by (14.11) . Thus we have shown the following result. 

Theorem 4.2. A DaS point process $ with the Levy measure supported by finite 
measures (equivalently, with a spectral measure o supported by Mi) can be represented 
as a cluster process with Poisson centre process on Mi driven by intensity measure a 
and with daughter processes being Sibuya processes Sib(a,/i), \i G M x . Its p.g.fl. is 
given by (14. 2p . 

As a by-product, we have established the following fact: since Sibuya processes are 
finite with probability 1, then the cluster process is finite or infinite depending on the 
finiteness of the Poisson processes of centres. 



G r [u] =E Y[ u{xi 

X;gT 
oo 

n=l 



1-fl 



u(xi) 



X r ' 
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Corollary 4.3. A DctS point process is finite if and only if its spectral measure a is 
finite and is supported by finite measures. 

If a = 1, then the Sibuya process consists of one point distributed according to 
fj, and the cluster process represents a Poisson process with intensity measure Jt(-) = 
f M fj,(-)a(dfi), which is clearly discrete 1-stable. 

By Corollary 13.61 £ is discrete a-stable if and only if £ can be represented as a sum 
of Poisson number of independent Sibuya distributed random variables, which is proved 
for discrete stable laws in [10]. 

Later we will make use of Theorem 14.21 for the case of an infinite countable phase 
space X. The finite case is considered below. 

Definition 4.4. Let fi = (/ii, , . . . , fid) be a probability distribution. A random vector T 
has multivariate Sibuya distribution with parameter measure \i and exponent a. 6 (0, 1] 
if its p.g.f. has the following form 

d d 

E J] zl- = 1 - (£(1 - z n )u. n y . (4.3) 

n=l n=l 

If d = 1 and /i — 1, then the multivariate Sibuya distribution becomes the ordinary 
Sibuya distribution with exponent a. For d > 2, the marginals T n of a multivariate 
Sibuya vector T = (T 1; , . . . , T d ) have the p.g.f. given by 

Ez^ = l-iu n (l-z n ) a . 

Thus T n takes value with probability 1 — fi n , but conditionally of being non-zero it is 
Sib(ct)-distributed. 

Example 4.5 (DaS random vectors). Let £ = ■■■,&) be a DaS vector. By Theo- 
rem [4721 it can be represented as a sum of multivariate Sibuya Sib(a, /Xj) vectors where 
fii are chosen from a finite Poisson point process on with some intensity measure a. 

Example 4.6 (Stationary DaS processes). As in Example 13. 9[ let X = M d with v being 
proportional to the Lebesgue measure and \x being the uniform distribution on the unit 
ball centred at the origin. Then the corresponding DaS process is a cluster process which 
can be described by the following procedure. First the Boolean model of unit balls in IR d 
is generated with the centres following the Poisson point process with intensity measure 
u, see [25] • Then a Sib(a) number of points is thrown into each such ball uniformly and 
independently from the other balls. The set of thus generated points is a realisation of 
the DaS process. 

In a more general model, the uniform distribution on the unit ball can be replaced 
by any probability distribution kernel P(dy,x), x G X. A typical realisation of such a 
process in R 2 with v being proportional to the Lebesgue measure on [0, l] 2 and P{dy, x) 
being the Gaussian distribution centred at x and with i.i.d. components with variance 
s 2 is presented in Figure [TJ The avoidance probabilities of the obtained point process 
are given by 

P{$(£?) = 0} = exp{- J (jMi(s-\B + x))J dx} , 
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Figure 1: A realisation of a DaS stationary process with the Gaussian kernel. Different 
shades of gray correspond to different clusters. 



where /xo is the standard Gaussian measure in M 2 . 

Figure [1] shows that the realisations of such a process are highly irregular. Many 
clusters with small number of points appear together with huge clusters, so that the 
resulting point process has infinite intensity measure. 



5 Regular and singular DaS processes 



We have shown in the proof of Theorem 13.31 that a DaS point process is necessarily 
infinitely divisible. It is known (see, e.g., 0, Th. 8.4.V]) that the p.g.fl. of an infinitely 
divisible point process admits the following representation 



G$[u] = exp 



M 



3 (logtt,¥>) 



1 



Q(d<p) 



(5.1) 
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where Ao is the space of locally finite non-empty point configurations on X and Q is a 
locally finite measure on it satisfying 

Q({<^ E Wo : <p(B) > 0}) < oo for all B E B {X) . (5.2) 

This measure Q is called the KLM-measure (canonical measure in the terminology of 
[T6] ) and it is uniquely defined for a given infinitely divisible point process $. Such 
point process $ is called regular if its KLM measure is supported by the set {tp 6 Mq : 
<p(X) < oo}, otherwise if is called singular. 

It is easy to see that the expression (13.71) for p.g.fl. of DaS process combined with 
(HOD has the form (ETTJ) with 

Q(-)= / P M (0A(dA0, (5-3) 
Jm\{o} 

where P M is the distribution of a Poisson point process driven by intensity measure /i. 
Moreover, the requirement (15.21) is exactly the property (12.71) of the Levy measure with 
h = l B . 

The following decomposition result is inherited from the corresponding decomposition 
known for infinitely divisible processes, see, e.g., [2J Th. 8. 4. VII]. Denote My = {// e 
M \ {0} : n(X) < oo} and M 00 = {/i6M: fi(X) = oo}. 

Theorem 5.1. A DaS point process with Levy measure A and spectral measure a can be 
represented as superposition of two independent infinitely divisible processes: a regular 
and a singular ones. The regular process is a cluster process having p.g.fl. given by (14. 2 j) 
with spectral measure <t|mi = ImJ and the singular process is a cluster process with 
p.g.fl. given by (\3.7\i corresponding to Levy measure AIm^ (or by (13. lip with spectral 
measure a|§\ M J- 

Proof. Represent A as the sum A| M/ + A| Moo of two orthogonal measures: A restricted 
on the set My and on M^. Then Q also decomposes into two orthogonal measures: 
the one concentrated on finite configurations and the one concentrated on infinite point 
configurations (because the corresponding Poisson process 11^ is finite or infinite with 
probability 1, correspondingly). Finally, G$[u] in (15 .ip separates into the product of p.g. 
fl.'s corresponding to the cluster process studied in the previous section and the second 
one with p.g.fl. (13. 7p with A replaced by A| Moc . □ 

Recall that a measure \i is called diffuse if u.({x}) = for any x G X. Representation 
( 15. 3p of the KLM measure immediately gives rise to the following results analogous to 
the properties of infinitely divisible processes. 

Theorem 5.2. A DaS point process $ with Levy measure A and spectral measure a is 

(i) simple if and only if A (respectively a) is supported by diffuse measures, c.f. /ITjl 
Prop. 2.2.9]; 
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(ii) independently scattered if and only if o is supported by the set {6 X : x G X} C § 
of Dirac measures, c.f. fffij Prop. 2.2.13]. 

(in) The distribution P x of a DaS processes $1 with spectral measure <r 2 is absolutely 
continuous with respect to the distribution P 2 of another DaS process $ 2 with 
spectral measure <r 2 and the same a if and only if there exists a measurable set 
A C Moo n S = §oo such that cxi | Soo = cr 2 | AnStx) and ai\ Ml < cr 2 | Ml - 

If X is compact, then all DaS processes are regular. In the setting of Example I3.9[ 
$ is singular if and only if [i is infinite. 

Consider the practically important case of a stationary DaS process $ in R d . Ac- 
cording to [5J Prop. 10. 3. XII], $ is mixing if and only if for all A, B G Bq 

Q{y? G No : y?(A) > and ip(B + x) > 0} — > as ||x|| — > 00 . 

In view of (15.31) chosing sufficiently large ||a;|| such that An (B + x) =0, this condition 
writes 

/ (1 - e^ (A) )(l - e -^ B+x) ) A(dfi) -> . 
Jm\{o} 

By (12. 7p and the dominated convergence theorem, $ is mixing if and only if fi(B+x) — > 
as 1 1 a; || — > 00 for all B G Bq and all fi from the support of the spectral measure a. This 
is clearly the case if a is supported by finite measures, i.e. for regular DaS processes. 
Similarly, the ergodicity (or weak mixing) condition from [5J Prop. 10. 3. XII] implies that 
$ is ergodic if and only if 

n~ d [ (1 - e-^ B+x) )dx -> as ||x|| -> 00 
for all B <E Bo and all /x from the support of a. 



6 Discrete stability for semigroups 

Let (S, ©) be an abelian semigroup with binary operation © and the neutral element 
e. We require that S possesses a basis X C S such that each element s £ S \ {e} is 
uniquely represented by a finite linear combination with positive integer coefficients 

s = mxi © • • • © n k x k , (6.1) 

where {xi, . . . , x^} Q X \ {e}. Here and below, ns stands for the sum s © s © • ■ ■ © s 
of n > identical elements s 6 S with convention 0s = e. The required uniqueness of 
representation (16 .ip implies that no element x of the basis can be expressed as a linear 
combination of others otherwise x would be multiply expressed as 1 • x, on one hand, 
and that other linear combination, on the other. It also rules out existence of cycles 
where the neutral element can also be obtained as a non-trivial sum of other elements. 

Note that {xi, . . . ,x k } with multiplicities {n 1; . . . ,n k } can be regarded as a finite 
counting measure on X, so that (16. ip establishes a map s 1— > <p between s G S and 
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a finite counting measure (p on X. Equip X with the topology that makes it a locally 
compact second countable space (if X is countable the conventional choice is the discrete 
topology). Define a er-algebra on S as the inverse image of the a-algebra on the space 
of counting measures under the map x i— > (p. A random ^-valued random element £ 
is defined with respect to the constructed a-algebra. Then the corresponding counting 
measure becomes a point process on X. 

The sum of random elements in S is defined directly by the ©-addition. The mul- 
tiplication t o £ of a random element £ by t G (0, 1] is defined as the random element 
that has the basis decomposition t o the latter obtained by independent thinning 
of $(£). By Theorem 13.11 and the uniqueness of the representation, the corresponding 
operation satisfies t\ o (t 2 o £) = (ii^) ° £ for £ij *2 G (0, 1] and 



for independent £ and Note also that t o e = e for all t G (0, 1]. 

The imposed distributivity property (16. 2\\ is essential to define the stability property 
on a semigroup (even for deterministic £ and £'), but it also rules out the possibility to 
extend the scaling operation to the cases where the decomposition (16. ip is not unique 
or the coefficients are allowed to be negative, which is the case if S is a group. Indeed, 
assume that a © b = e for non-trivial a,b G S. Then t o (a © b) = e, while t o a © i o 6 
is the sum of two non-trivial independent random elements and so cannot be e almost 
surely. This observation explains why the classical notion of discrete stability cannot be 
extended to the class of all integer valued random variables. 

Similar difficulties arise when defining discrete stability for the max-scheme being an 
example of an idempotent semigroup. Indeed, if S is the family of non-negative integers 
with maximum as the semigroup operation, then a © a = a. After scaling by t G (0, 1) 
we obtain that t o a is the maximum of two independent random elements distributed 
as t o a, which is impossible. 

Definition 6.1. An S- valued random element £ is called discrete stable with exponent 
a (notation: DaS) if 



for any t G [0, 1], where £" and £ are independent identically distributed. 

Thus, £ is DaS if and only if the corresponding point process $(£) is DaS, see 
Definition 13.21 Thus, the distributions of DaS random elements in S are characterised 
by Theorem 13.31 Since the basis decomposition (16.11) is assumed to be finite, $(£) is 
a finite point process, meaning that the spectral measure a is finite and supported by 
finite measures, see Corollary 14.31 

To formalise the presentation, consider characters \ on S, which are homomorphisms 
between S and [—1,1] with the multiplication operation. While in general characters 
are homomorphisms between a semigroup and the unit complex disk, the semigroup S 
is equipped with the identical involution, implying that characters are real-valued. It is 




(6.2) 




(6.3) 
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shown in [7j that if the characters form a separating family S, then the distribution of 
£ is uniquely defined by the expectations E%(£) for x G S 1 . It is easy to see that 

E x (0=G m [ X }- 

Note that 

E X (t°0 = G* (0 [l-t(l-x)]. 
Corollary 13.41 yields that £ is DaS if and only if 

E X (0 = exp |- jf <1 - *, //> } , 

where § is defined in (12.41) . 

A random element v in S is said to have Sibuya distribution if $(u) is a Sibuya point 
process on X. Theorem 14.21 immediately implies the following result. 

Theorem 6.2. A random element £ in S is DaS if and only if £ is can be represented 
as the sum of a Poisson number of i.i.d. Sibuya random elements. 

The LePage representation of DaS processes from Corollary 13.51 implies that £ is DaS 
if and only if it can be represented as the sum ^i>i ^Hk wnere £i, £2, • • • are i-i-d. 

and the sum is almost surely finite. 

The case of a finite X is important for semigroups with a finite basis. Then (16. ip 
establishes a homomorphism between S and Zl, i.e. the family of (i-dimensional vectors 
with non-negative integer components and addition as the semigroup operation. Thus, 
all random elements in semigroups with finite basis can be treated as random vectors 
in Z+. Its discrete stability property has been addressed in Example 14.51 

Now consider the case of a non-finite countable X considered a basis for a discrete 
semigroup. 

Example 6.3 (Arithmetical semigroup). Consider an arithmetical semigroup which is by 
definition a commutative semigroup G generated by a countable subset V = {pi, P2, • • • } 
called the generalised primes such that any element g G G admits a unique factorisation 

9 = \{v np , 

pev 

where n p are non-negative integers only finite number of which are non-zero. The sum 
^2 n p defines the degree d(g) of g and it is assumed that the total number of elements g 
with d(g) < x is finite for all x > 0, see, e.g., [T5]. Examples include the set of positive 
integers with the prime numbers as the generator set V, or, more generally, Beurling's 
generalised prime numbers with the so-called Delone property which uniformly bounds 
(from above and from below) the gaps between the successive elements of G C M+ and 
implies uniqueness of the factorisation, see [2], [8]. 

To any finite point process $ on V there corresponds a positive G-valued random 
variable 

per 
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and the nitration operation naturally extends to it: t o is, by definition, £ to $, so that 
t o £ is obtained by independently taking out the generalised prime divisors of £ with 
probability 1 — t. Consequently, a class of discrete multiplicatively stable distributions 
can be defined as the distributions of G-valued random variables £ satisfying 

(t 1/a oe')-((i-o 1/a °e") = e, 

where, as before, £" are independent copies of £. Then £ = is multiplicatively 
a-stable if the corresponding process $ on V has p.g.fl. given by (13. lip . 

Assume now that V C [1, oo) as it is true for prime numbers (and also for Beurling's 
generalised primes). Extend the domain of the p.g.fl. to the class of monotone pointwise 
limits of the functions u such that 1 — u G BM(V) allowing for infinite values of the 
p.g.fl. and consider h s : V \— > (0, 1) such that h(p) = p s for some s. Then 

E& = ETTp'*W = G 9 [h.] = exp{- / (1 - V P > P )V(^)) , 

where Mi is the set of probability distributions on V. The expression above is finite at 
least for all s < and can be used to numerically evaluate the distribution of via the 
inverse Mellin transform. 

Now consider the case of independently scattered $, where a is concentrated on de- 
generated distributions (see Theorem I5.2l( ii)) and thus can be identified with a sequence 
{cr p , p G V}. Then $ is a sequence indexed by p G V of doubly stochastic Poisson 
random variables v^ v with parameters ( p which are independent positive StaS random 
variables with Laplace transform E e" ftp ^ = exp{— hp<j p }. Now the distribution of a 
can be explicitly characterised: if n = Yl p evP kp j then 

= n} = P{ % = k p Vp G V} = J] E [^e-<"" . 

per p ' 

In particular, for any q G T 3 we have that 

P{£* = <?} = II P K = °> E [C* e_C '] = e~° iv)+ ^ E[C ? e"^] . 

Take, for example, a = 1/2. Then the density of ( p is given by 

/cp(t) = vS^ exp{_ap/(4t)} ' £ -° 

leading to 

Example 6.4. Consider the family S of all finite abelian groups with the semigroup 
operation being the direct product. The main theorem on abelian groups states that 
each such group can be uniquely decomposed into the direct product of cyclic groups 
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with orders being prime numbers and their natural powers, see, e.g., [TT1 Th. 7.2]. 
Thus, the basis X is the family of prime numbers and all powers of prime numbers. 
The multiplication of a cyclic group of order p by real number t 6 (0, 1] is defined as a 
random group where each factor from its decomposition is eliminated with probability 
1 — t. A spectral measure defined on the sphere in X then determines the distribution 
of a random stable finite group. 
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